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ABSTRACT 

We investigate the transition to perturbative QCD in Compton processes, by concen- 
trating on the specific reactions 717 — ► tc-/ and pion photoproduction at moderate energy 
scales. New sum rules for each of the helicities involved in the scattering and the cross sec- 
tion are given, together with a detailed stability analysis. Our results are then compared 
with those derived from higher-power factorization incorporating Sudakov suppression 
for quark elastic scattering. An overlap of predictions obtained by the two approaches 
is observed at a momentum transfer as low as 4 GeV 2 and at a photon scattering angle 
around 40°. Our work shows that factorization theorems and sum rule methods com- 
plement each other in the description of Compton scattering at moderate energy scales, 
and give information on the transition to perturbative QCD in these processes. The sum 
rule formalism is further applied to the crossed reaction 77 — > 7r + 7r~, and is in very good 
agreement the experimental data. 
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1 Introduction 



Exclusive processes have been one of the most challenging testing grounds for perturbative 
QCD (PQCD). Although the description of exclusive processes using PQCD is considered 
successful for momentum transfer Q 2 going to infinity, it remains a controversy whether 
PQCD is applicable to exclusive processes at moderate Q 2 . This issue has been widely 
discussed in the literature, however, general agreement on it is still missing. It has been 
shown that the standard leading-order and leading-twist factorization [|I|] fails to give 
reliable predictions to hadron form factors at experimentally accessible Q 2 because of the 
dominance of soft contributions from the end-points of parton momentum fractions [0. 
More recently, the understanding of exclusive processes has been improved by studying 
the transverse momentum dependence of scattering amplitudes || , which is introduced by 
transverse momenta carried by the collinear partons that enter hard scatterings. Modified 
perturbative expressions for pion and proton electromagnetic form factors including this 
higher-power dependence have been given [§], . It is found that the all-order summation 
of the higher-power effects leads to suppression of end-point contributions, and extends 
the applicability of PQCD down to 2-3 GeV. That is, hadron form factors transit from 
non-perturbative to perturbative QCD around this scale. 

However, it is still unclear how far low in Q 2 the modified perturbative picture is yet 
successful for other more complicated exclusive processes. In this paper we shall inves- 
tigate the transition to PQCD in Compton processes by comparing predictions from the 
above modified perturbative formalism with those from QCD sum rules @, a method to 
study non-perturbative properties of hadrons. A sum rule approach to a simple Compton- 
type process, pion Compton scattering, has been discussed in refs. and 0, in which 
the similarity between fixed-angle Compton scattering and form factors was explored. 
Sum rules for the sum H of the two helicity invariant amplitudes, Hi and H 2 , involved 
in pion Compton scattering was given in ||, and preliminarily compared to the modified 
perturbative predictions in ||. A transition scale at Q 2 as low as 4 GeV 2 , for photon 
scattering angle around 40°, has been observed. Compared to the sum rule formalism 
for form factors, the analyticity region of Compon scattering is more severely constrained 
[0, §, 10 j. In order to maintain the stability of the sum rule for H in the finite analyticity 



region, we have proposed a modified phenomenological parametrization for resonance and 
continuum states in [p~T|l . 

PQCD and sum rule approaches employ very different physical pictures: the former 
takes a hard-scattering view of exclusive processes, while the latter is described by the 
dominance of "Feynman mechanism" fT2 |. It has been pointed out |4], [7]] that the PQCD 
approach should be complemented by sum rule methods in order to fully understand the 
behaviors of simple elastic scatterings at moderate energy scales. The parallel discussion 
based on both PQCD and sum rules is of interest, since it can help to better understand 
the transition to PQCD in exclusive processes, and clarify, to some extent, the questioned 
applicability of PQCD. 
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In this paper we shall investigate individual sum rules, for Hi and H 2 respectively, 
in great detail, together with their modified stability analysis ||1 1|| . The cross section of 
pion Compton scattering is also given. All the steps are then repeated using the modified 
PQCD formalism, and predictions from the two approaches are compared. We find that 
perturbative contributions dominate at large angles, while at smaller angles sum rule 
contributions are important. The overlap between the two descriptions appears at about 
40° of photon scattering angle and at Q 2 around 4 GeV 2 , similar to the behavior observed 
in the sum rule for H [[|. We then confirm the conclusion drawn in |0| that large-angle 
pion Compton scattering can be treated by pQCD reliably. To justify our formalism, the 
sum rule methods are further applied to the crossed version of Compton scattering, pion 
photoproduction 77 — > 7r + 7r~, for which experimental data are available It is found 
that our predictions are consistent with the data. 

This paper is organized as follows. A brief illustration of the sum rule formalism is 
given in section 2, where the main features of the sum rule for H, extensively studied in 
H [|, are reviewed. We present in section 3 the complete evaluation of the individual 
sum rules for Hi and H 2 , which include both the local duality contribution and the 
power corrections. The stability analysis of the new sum rules based on the modified 
phenomenological model [11] is performed. H and H\ are evaluated using the PQCD 
formalism in section 4. Section 5 contains numerical results of the two helicities and the 
cross section derived by the two approaches, and we specify the transition region. The 
sum rule for pion photoproduction is also analyzed, and the results are compared with the 
data. Section 6 is the conclusion. Three appendices are inserted to clarify the conventions 
and to decribe some details of the calculation. 



2 Sum Rules for H 

As the relevant energy scales approach a resonance region, a new complexity in the de- 
scription of exclusive processes shows up, because non-perturbative effects, related to the 
dynamical behaviors of QCD vacuum, can no longer be neglected. Such effects can not be 
incorporated into a direct perturbative treatment based on factorization theorems, and 
to keep them into account one must resort to completely different methods, the most 
successful one being QCD sum rules. Sum rule methods have been recently extended 
to fixed-angle Compton scattering from their usual application to two- and three-point 
processes in 0. So far, all previous works have concentrated on the specific combination 
of the invariant amplitudes, H = Hi + H 2 . We have communicated before in a brief letter 
H that sum rule predictions for H dominate at intermediate photon scattering angles 
over perturbative calculations. 

QCD sum rules connect the timelike region of a suitably chosen correlator to the space- 
like region in the complex planes of external virtualities p 2 by a dispersion relation. For 
pion Compton scattering the relevant correlator corresponds to the lowest-order diagrams 
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without virtual gluons, as shown in fig. 1. In the case of four-point functions the existence 
of an analiticity region for the correlator in the pf planes, at fixed Mandelstam invariants 
s, t and u, is far from being obvious. 

In the deep Eucliden region, the spacelike region with large negative virtualities, the 
operator product expansion (OPE) for the correlator is usually valid, and can be cal- 
culated using perturbation theory. The result is organized in terms of a lowest-order 
perturbative contribution plus power corrections, the latter being parametrized by the 
lowest dimensional QCD condensates. These corrections can be obtained (to lowest or- 
der in a s ) in two standard ways: by Cutkosky rules, and by Borel transforms |7J. Both 
methods (TB| |T2j have been extensively employed for three-point functions. The regions of 
analyticity in the case of massive and massless correlators have been described in detail in 
ref. [TUj. Radiative corrections to the dispersion relations, which enter in the description 
of these processes, are also calculable by the methods discussed in ref. [|K] . 



The timelike region of the correlator cannot be described by perturbation theory for 
the virtualities close to the lowest resonances, and the residues at the resonance poles 
are known only up to their symmetry properties. Therefore, a phenomenological model 
parametrizing the resonance and the continuum contributions has to be proposed. The 
continuum contribution is usually chosen as the perturbative part on the Operator Prod- 
uct Expansion side of the sum rules. 

For forward or backward scattering the Compton processes cannot be described by 
the extended sum rule formalism. A simple observation of this fact is that the involved 
leading spectral function becomes singular as the Mandelstam invariants s and t vanish 
[|i~0f . Hence, the OPE is applicable only in a suitably chosen angular range. In order 
not to include extra w-channel singularities, we study the correlator in a finite region in 
the pf planes. Even for a finite analiticity region, the existence of a dispersion relation 
guarantees that the behavior in the timelike region can be related to that in the spacelike 
region. This relies on the introduction of a modified Borel transform , with finite radius 
A 2 , to characterize the analyticity region. A spurious dependence on the "Borel radius" 
A 2 is then brought into the dispersive representation of the coefficients of such corrections. 
As discussed in jTl|, this A-dependence is unphysical and spoils the stability of the sum 



rule. 

The coefficients of the OPE are usual Feynman integrals, now in a dispersive form 
||. It is possible, however, to extend the dispersive representation of each Feynman 
integral (and also of the power corrections) to all positive values of Si and s 2 , in the 
case of massless correlators. For pion photoproduction the pion pole is approximately 
set at p\ = p\ = 0. If one wants to extract the contribution of additional states from 
the phenomenological side of the sum rule or, for instance, extend the method suggested 
in ref. Jfj to the proton case, then the spectral functions evaluated by the OPE have 
an explicit mass dependence and additional singularities compared to the massless case. 
These singularities, specifically, are described by Landau surfaces and disappear as the 



mass in the correlator is set to be zero. The discussions of all these issues is in ref. ITU 
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From a physical viewpoint, it is expected that at large virtualities the spectral densities 
should give negligible contribution to the dispersion integral. It is also easy to show that 
most of the contribution to the spectral density comes from the region of approximately 
equal virtualities s\ and s 2 (si ~ s 2 ). 

Based on these arguments, we have proposed a modified phenomenological model for 
resonance and continuum states in [TT| to avoid the undesirable A 2 dependence. This new 
ansazte removes the contributions of large virtualities both from the phenomenological 
side and from the OPE side of the sum rules. In other words, according to this model for 
the continuum (on what is termed "the phenomenological side" of the sum rule), all the 
coefficients of the OPE cannot contribute to the sum rule for large virtualities. 



Application of this ansazte to the pion form factor |TTJ gives results which are not 



substantially different from what have been obtained in the literature [[12], |I5 |. A stability 



analysis of the sum rule for H based on the modified phenomenological model has been 



performed in |TT[, where the local duality interval Sq and the Borel mass M are found 
to take values very close to those in the form factor case [|T^, [T5| . In order to make our 
discussion self-contained, we briefly summarize the derivation of the sum rule for H 
We then compute H\ and H 2 , from which the cross section is derived accordingly. 
We start with the following four-point correlator 

^a^x{p\,pl, s,t) = i J d A xd A yd 4 zexp(-ip 1 ■ x + ip 2 ■ y - iqi ■ z) 

x(0|T( ?7(J (2 / )J M ^)J l/ (0)4(x))|0} ) (1) 

where 

2 1 - 

Jn = 2^7m m - g^ci, r] a = vr/5lad (2) 

are the electromagnetic and axial currents, respectively, of up and down quarks. The 
on-shell momenta q\ and q 2 are carried by the two physically polarized photons. The two 
pion momenta are denoted by p\ and p 2 , with s\ = p\ and s 2 = p\ the virtualities. The 
Mandelstam invariants s = (pi + q\) 2 , t = (p 2 — Pi) 2 and u = (p 2 — q\) 2 obey the relation 
s + t + u = S\ + s 2 . The invariant amplitudes Hi and H 2 involved in pion Compton 
scattering are then obtained through the expansion of the matrix element 



Mfa, = i J d A ze-^ z (p 2 \T (Jp(z)J v (0)) |pi) (3) 
for a specific time-ordering as 

ilV = HiMeWeV + H 2 (s,t)eW e V . (4) 

The helicity vectors e'- 1 '' and are defined in [[?], [TTH , satisfying the orthogonality con- 
dition ■ = —Sij. The expression of the sum rule for H is derived by contracting 
Tcr^A with —g^ u n a n x , where 

(e®±ieWy (5) 
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is a suitable projector in analogy to that for pion form factor |L2|. The phenomenological 
model employed in sum rules are characterized by a resonance contribution from the 
double poles of the pion states, and by a continuum one A pert |7], |j for pf > s , which is 
exactly the perturbative contribution on the OPE side: 

A (p 2 i,s,t) = f%n-p l n-p 2 (2ir) 2 5(pl-ml)5(pl-ml) 

xH(s,t)+A^ipls,t)[l-9(s Q -pl)9(s -pl)) . (6) 

A modified Borel transform [[F] 

B' = j M j ^% e - P l/M? e -pl/Ml ^ _ e -(^- P i)/M^ ^ _ e -{\*-pl)/M^ (7) 

was then proposed to define the analyticity region, with C a contour of finite radius A 2 , 
which is allowed to vary from sq to (s + t)/2. The factor 1 — exp[— (A 2 — pf)/M 2 ] in 
(|7|), which is new compared to the standard Borel transform |12], [15], gives the constraint 
pf < A 2 , and thus excludes the w-channel resonances from the spectral density, ensuring 
that the transform is still regular when C crosses the branch cuts in the pf planes. This 
modification introduces the extra unphysical parameter A 2 , in addition to the usual Borel 
mass M 2 . Hence, the stability of the sum rule has to be established in the two-dimensional 
M 2 -A 2 plane. Using eq. (|6|), the sum rule for H has been derived [0], in which the 
argument Si of the perturbative spectral density on the OPE side runs up to so, due to 
the cancellation by the continuum contribution on the phenomenological side; while 
in the power corrections runs up to A 2 . A naive numerical analysis based on the above 
formalism has been performed, and results show a strong sensitivity to the variation of 



A 2 JIJ. 



The strong A dependence of the sum rule is not desired, since A is an unphysical 
parameter, and a physical quantity like the invariant amplitude should be insensitive 
to it. Therefore, we took an alternative approach to avoid the ^-channel singularities. 
We have proposed a modified phenomenological model with the continuum contribution 
A pert replaced by A OPE , the full spectral expression on the OPE side of the sum rule. 
This modification makes sense, because the region with large virtualities pf > sq can 
be regarded as perturbative, and an OPE is allowed. With this choice we are requiring 
that the phenomenological parametrization for the continuum from s to A 2 truncates not 
only the perturbative part, but also the power corrections, on the OPE side. Since the 
contributions from the large-virtuality region (sq, A 2 ) have been removed by the above 
cancellation, s^'s never reach the upper bound A 2 , and the remaining A dependent factors 
1 — exp[— (A 2 — pf)/M 2 ] from the transform (|7]) can be dropped. Therefore, the overall 
dependence on the Borel radius disappears completely from the sum rule. In conclusion, 
we suggested a modified phenomenological model with the original Borel transform 

, / . ^ pV Mt e - pV Ml (g) 

Jo M\ Jc Ml v ' 
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instead of a new transform with the original model as in JF). The resulting asymptotic 
expression for H based on the modified model is written as 



,,, ^ / s(s + t) 



U 2 H(s,t) 



t 

' r d S1 r ds 2P pcvt + —(g 2 ) r d Sl r ds 2P siu ™ 

\Jo JO 71 JO JO / 

+C quark vra s ((^) 2 ) (9) 
for large invariant Q 2 H], 



e -{s 1 +s 2 )/M 2 



Q 2 =~(s 1 + s 2 -t+ ^ (s 1 + s 2 - t) 2 - 4: Sl s 2 ^ , (10) 

where the perturbative, gluonic and quark contributions are given by, respectively, 

2560Q 14 t^(s,Q 2 , Si ,s 2 ) 
P " ^{s-2Q 2 ){AQ i -s 1 s 2 f{s 1 s 2 -2QH) ' 

20480Q 22 (s - 2g 2 )rS luon (g, Q 2 , s u s 2 ) 

27s(2Q 2 - Sl ) 2 (2Q 2 - s 2 ) 2 (4Q 4 - Sl s 2 f{2Q 2 s - sis 2 )(4Q 4 - 2QH + Sl s 2 ) 2 ' 
rk _ 16 (8MV + 8M 4 st + 8M 2 s 2 t + 8M 2 st 2 + AsH 2 + Ast 3 + t 4 ) 



with 



T P ert = ( s -Qiy(2Q A ss 1 -Q 2 s 2 s 1 + 2Qhs 2 -Q 2 s 2 s 2 
—2Q 4 s 1 s 2 - QQ 2 ssis 2 + 3s 2 s 1 s 2 ) , 

r gluon = _ g Ql2 s _ gglO s 2 + ggg8 s 3 _ g 4 g6 s 4 + 16 g4 s 5 

+8Q 10 ssi + 8Q 8 s 2 si - 108Q 6 s 3 si + 104QVsi - 26QVsi 
+8Q 10 ss 2 + 8Q 8 s 2 s 2 - W8Q 6 s 3 s 2 + 104Q 4 s 4 s 2 - 2QQ 2 s 5 s 2 . 



'12) 



Approximate methods for the evaluation of power corrections have been developed in 
ref. ||. A further analysis of the diagrammatic expansion jnj shows that it is possible 
to take the limit of A — > oo in the dispersive representation of the coefficients of the 
power corrections, since their the spectral densities are globally well defined for massless 
correlators. The rigorous proof of these statements is discussed in ref. ||10|| . Note the 
upper bound s instead of A 2 in the integral for the gluonic power correction due to the 
cancellation from the phenomenological side. The gluon and quark condensates, (G 2 ) and 
((tjjijj) 2 ) , take the values 

— (G 2 ) = 1.2 x 10- 2 GeV 4 

71 

a s {{i)i)) 2 ) = 1.8 x KT 4 GeV 6 . (13) 
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Since the A dependence is removed completely, the stability analysis of (§) is straight- 
forward following a method similar to |12], |T5| . We concentrated simply on the variation 
of H with respect to M 2 . It has been found that as So = 0.6 GeV 2 there is the largest 
M 2 interval, in which H is approximately constant. Therefore, Sq = 0.6 GeV 2 is the best 
choice which makes both sides of the sum rule most coincident. This value of the duality 
interval is close to that given in the form factor case, and consistent with its conjectured 
value of 0.7 GeV 2 in |§. Different sets of s and t have been investigated. The best value 
of So does not vary significantly, and H is almost constant within the range 2 < M 2 < 6 
GeV 2 . 

Results for H at different photon scattering angles 9* in the Breit frame, sin(#*/2) = 
—t/ (s — u), with so = 0.6 and M 2 = 4 GeV 2 have been obtained ||, |H|| . Basically, they 



show a similar dependence on angles and momentum transfers \t\ to those derived using 
local duality approximation |5| . These predictions have been compared to the perturbative 
predictions obtained from the modified factorization formula 0, which will be discussed 
in detail below. Sum rule results are always larger than the perturbative results at smaller 
angles. The transition to PQCD appears at about \t\ = 4 GeV 2 and at 9* = 40°, where 
the perturbative contributions begin to dominate. 



3 Sum Rules for H\ and H2 

Following the formalism outlined in section 2, we derive sum rules for the helicity invariant 
amplitudes H\ and if 2, which are extracted by contracting T a ^ v \ with e^^e^ u n a n x and 
e (2)M e (2)f n °" n A respectively. Similarly, the modified phenomenological model for Hi is given 
by 

Ai(p 2 ,s,t) = f 2 n-p l n-p 2 (2ir) 2 5(p 2 1 -ml)5(pl-m 2 7T ) 

xH l (s,t) + Af PE (p 2 ,s,t)[l-9(s -p 2 1 )9(s -pl)] . (14) 

The method to evaluate the perturbative spectral density and the power corrections from 
quark and gluon condensates is described in M. The result can be organized in the form 



//#i(M) 



's(s+ty 



-t 



\J0 JO 7T JO JO J 

+cr vk ™ s ((^) 2 ) . (15) 

The perturbative and gluonic contributions for Hi are given respectively by 



-1280Q u rr\s,Q 2 ,Si,s 2 ) 
Pl 3 7T 2 (s-2g 2 )(4Q 4 -s 1 s 2 ) 5 (s 1 s 2 -2g 2 s) 
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gluon 



Pi 



with 



20480Q 22 (s - 2Q 2 )rf luon ( S , Q 2 , s u s 2 ) 



27s(2Q 2 - Sl ) 2 (2Q 2 - s 2 ) 2 (4Q 4 - Sl s 2 f{2Q 2 s - Sl s 2 )(4Q 4 - 2Q 2 s + Sl s 2 ) 2 ' 

(16) 



r P-t = (s _ g2 )2(2 g4 SSi _ Q 2 s 2 Si + 2Q 4 S$2 _ 

-(2Q 6 s - 13QV + 12QV - 3s 4 )sis 2 , 

r gluon = _ 36 QlO s 2 _ U g8 s 3 + 32QV - 8QV 

+22Q 10 ssi + 2Q 8 s 2 si + 99Q 6 s 3 si - lOOQVsi + 25QVsi 



+22Q w ss 2 + 2QVs 2 + 99QVs 2 - 100QVs 2 + 25QVs 2 



(17) 



The quark power correction for Hi is computed as 



Arrk ^ 256Q 6 (2Q 2 - g )rr ' k 



(2Q 2 - pi 2 ) 2 pi 6 (2Q 2 - p 2 2 ) V (4Q 4 - 2Q 2 s + Pl 2 p 2 2 ) 



where 



.quark = ^ g6 _ gg4 s + ^2 _ ^4^2 + ^2^2 _ 3^2 _ ^4^5 

+6Q 2 sp 2 2 - 3s V) (pi 4 +p 2 4 



(19) 



From eq. fll8"|) £7^ uark i s obtained through the application of the standard Borel transform 
in eq. (||) as 



Cr rk = B{p\^Ml)B{pl^Ml)Nr k \ MuM2=M 

16 (8M 4 s 2 + 8M 4 st + 16M 2 s 2 t + 16M 2 st 2 + 4s 2 t 2 + 4st 3 + t A ) 
~~9 MH 3 

The corresponding quantities associated with H 2 are given by 
1280Q 14 r 2 pcrt ( g ,Q 2 , gl , g2 ) 

P2 - 



• (20) 



gluon 
p\ 



3tt 2 (s - 2Q 2 )(4Q 4 - Sl s 2 f( Sl s 2 - 2Q 2 s) ' 

20480Q 22 (s - 2Q 2 )rf luon ( S , Q 2 , s u s 2 ) 



27s(2Q 2 - Sl ) 2 (2Q 2 - s 2 ) 2 (4Q 4 - Sl s 2 fi2Q 2 s - Sl s 2 )(4Q 4 - 2Q 2 s + Sl s 2 ) 2 ' 

^yquark _ 128 s(.S + t) 

° 2 " ~9~ltf¥~' 
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with 



r P ert = 3(s-Q 2 ) 2 (2Q 4 ss 1 -Q 2 s 2 s 1 + 2Q 4 ss 2 -Q 2 s 2 s 2 ) 
-(4Q 8 + 6Q 6 s - 39QV + 36QV - 9s*) Sl s 2 , 

r gluon = _ 8Q 12 S + 28 glO s 2 + g2 g8 s 3 _ gQ Q 6 ^ + ^qA ^ 

-14Q 10 ssi + QQ 8 s 2 Sl - 207Q 6 s 3 si + 204Q 4 s 4 Si - 51QVsi 
-14Q 10 ss 2 + QQ 8 s 2 s 2 - 207Q 6 s 3 s 2 + 204Q 4 s 4 s 2 - 51QVs 2 . 



[22) 



It can be easily checked that the combination of the above two sum rules, Hi + H 2 , is 
equal to that of H as shown in section 2. 



A similar stability analysis of eq. (15) is performed, and the variation of Hi with 
respect to M 2 for s = 0.5-0.7 GeV 2 at s = 20 and \t\ = 10 GeV 2 is shown in fig. 2. The 
region on the right-hand side of vertical bars is the one where the power corrections do 
not exceed 50% of the perturbative contribution. Obviously, in both of fig. 2a and 2b the 
curve corresponding to Sq = 0.6 GeV 2 exhibits a largest M 2 interval, in which Hi remains 
almost constant. By varying s and t, we find that the best choice for the duality interval 
is sq = 0.6 GeV 2 and the corresponding stability region is 2 < M 2 < 6 GeV 2 , the same 



as those for H JTT]]. We shall evaluate Hi by substituting s = 0.6 and M 2 = 4 GeV 2 into 
eq. ( |T5| ) in section 5 in order to compare it with predictions derived from the modified 
perturbative expression. The detailed calculation of all the relevant quantities in the sum 
rule approach is given in Appendix A. The full expressions for the perturbative spectral 
functions of Hi and H 2 are exhibited in Appendix B. 



4 The PQCD approach 

The modified pQCD formalism including transverse momentum dependence of the hard 
scattering, which will be presented below, has been discussed extensively in refs. ||, |5|, |§ . 
We shall summarize the basic reasoning leading to the modified perurbative expression, 
and review the calculation of the invariant amplitude H briefly [0. We then derive the 
first helicity amplitude Hi using the same method. The second invariant amplitude is 
obtained by H 2 = H — Hi. 

We consider the factorization formula for Hi in the conjugate configuration, or b, 
space due to the introduction of transverse momenta, b being the separation between 
two valence quarks. The perturbative expression includes the infinite summation of the 
higher-power effects at lower momentum transfers, which suppresses non-perturbative 
contributions. This effect, termed "Sudakov suppression" , exhibits a falloff of the invariant 
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amplitudes at fixed b faster than any power of Q. The modified perturbative expression 
reduces to the standard one as predicted by the dimensional counting rules [1J at very 
large Q 2 . The higher-power effects appear as large logarithms resulting from radiative 
corrections to the process, which will be analyzed below. 

Basic diagrams for pion Compton scattering in the PQCD approach are shown in fig. 3, 
which differ from fig. 1 by an extra exchanged gluon. In most cases, large logarithms from 
radiative corrections in exclusive processes do not cancel, and have to be summed up. 
The radiative corrections to fig. 3 can be classified into two categories: reducible and 
irreducible. The former, with the extra gluon connecting two incoming (or outgoing) 
valence quark lines, have both collinear and soft divergences. The latter, with the extra 
gluon connecting the incoming and outgoing quark lines, are lack of collinear divergences. 
Therefore, only the reducible radiative corrections are leading and should be resummed 
into the Sudakov factor mentioned above. If the main contributions are due to soft gluon 
exchange, the case in which the running coupling constant a s becomes large, lowest-order 
PQCD will not be reliable. With the introduction of the extra b dependence, a s remains 
small as long as b is small. The effect of Sudakov suppression provides the exact mechanism 
of confining the scattering process in the short-distance region. Therefore, the modified 
PQCD formalism is relatively self-consistent. The method to calculate pion Compton 
scattering based on these diagrams is similar to that developed for electromagnetic form 
factors [§. 

We start with the factorization formula for H(s,t), keeping transverse momentua 
carried by valence quarks in the pion wave function if) and the hard-scattering kernel Th, 

2 ,1 

H(s,t) = V / dx 1 dx 2 d 2 k Tl d 2 k T2 ip(x 2 ,k T2 ,p 2 ) 

xT Hl (xi,s,t,k T .) ip(x u k Tl ,pi) . (23) 

If kx in Th is assumed to give higher-power (k^/Q 2 ) correction and thus neglected, the k^ 
integration can be performed, and eq. (|23"D leads to the standard factorization formula [|TJ . 
However, this approximation is not proper when the exchanged gluon becomes soft. The 
contribution to the hard scattering from each diagram in fig. 3, obtained by contracting 
the two photon vertices with — is listed in table 1. All the contributions can be 
grouped into two terms (1 = 1, 2) using the permutative symmetry. 

Rewriting eq. (|23| ) in terms of the Fourier transformed functions, and inserting the 
large-6 asymptotic behavior of the wave function ||, we have 

Jjj rl roc 

H(s,t) = y~] / d^idx2 <p(xi)4>(x2) I bdbTnAxi, s,t,b) 
l=1 Jo Jo 

xexp[-S(x h b,Q)]. (24) 

where b, introduced by the Fourier transform, is the separation between the two valence 
quark lines as stated before. In the PQCD approach the pions are assumed to be on-shell, 
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pf = 0. Hence, we have t = —2Q 2 , as obtained from eq. (|10"D by setting Si to zero. Note 
the extra Sudakov factor exp(— S) compared to the standard factorization formula, which 
arises from the all-order summation of the collinear enhancements in radiative corrections 
to fig. 3. The exponent S is written as 0] 



2 2 w 

S(x 1 ,x 2 , b, Q, w) = J2( s ( x i' & > Q) + s ( 1 ~ x i' h iQ))~ ^~ ln — - 

Pi —o 



(25) 



i=i 



with 



s(£,b,Q) 



A (2) 



m \-b 



2/3i 



(q + b) 



A® fa 



-q 



ln(-26) + 1 ln(2g) + 1 



-b 



16/3? 

A® A® le 



27-1 ' 



In 



-b. 



' 32/3? 



ln 2 (2g) - In 2 (-26) 



(26) 



The variables q, b and w are defined by 



q = In [^/A] 
b = ln(6A) 
w = ln(w/A) , 



(27) 



where the scale parameter A = Aqcd will be set to 0.1 GeV. The coefficients /3j and A® 
are 



33 -2n 



f 
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153-19n/ 



24 



67 7T 
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A (D = 2 ; A (2) = _ _ _ _ 21 n + _^ l n _ 



9 



27 



:28) 



with n/ = 3 the number of quark flavors and 7 the Euler constant. The Sudakov factor 
is always less than 1 as explained in 0J, and decreases quickly in the large-6 region. The 
function 0, obtained by factoring the Q and b dependences from the transformed wave 
function into Sudakov logarithms, is taken as the Chernyak and Zhitnitsky model [pl| 

15/, 



<P CZ {x) 



v/2AT 



x(l-x)(l-2xY 



(29) 



where N r = 3 is the number of colors. 
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The transformed hard scatterings Tjy are given by 



~ 16nC F (e 2 u + e 2 d )a s (wi) 

= ~ "~ TT^ ~ K 



(1- Xl )(l-X 2 ) 

[(1 - Xl )t + s][(l -x 2 )t + s 
a 2 




[(1 -xi)t + u][(l -x 2 )t + u] 



4 l-x 2 



(30) 



from the classes of fig. 3a-3c, and 

T H2 = 32irC F e u e d a s (w 2 ) 



6(-r 2 )K [yj\r 2 \bj + e(r 2 ) l ^H^ (0^6) 

(1 + X 2 — XiX 2 )t 2 + (1 + x 2 — Xi)ut 

x 2 (l — Xi)s 2 
(1 + X\ — X\X 2 )t 2 + (1 + X\ — x 2 )st\ 
xi(l-x 2 )u 2 J 



X 




x 2 (l - x 2 ) 
from the classes of fig. 3d-3e with 



ri = xix 2 t, r 2 = x\x 2 t + x\u + a^s . 



(31) 



(32) 



Kq and Hq 1 ^ in eqs. (^) and ( |3~TD are the Bessel functions in the standard notation. The 
imaginary contribution comes from the case in which the exchanged gluons in fig. 3d and 
3e are on-shell, or r 2 vanishes. The argument wi of a s is defined by the largest mass scale 
in the hard scattering, 

i,t), w 2 = max (yj\r2\ , ^ 

As long as b is small, soft r\ does not lead to large a s . Therefore, the non-perturbative 
region in the modified factorization is characterized by large b, where Sudakov suppression 
is strong. Eq. (pl[) , as a perturbative expression, is thus relatively self-consistent com- 
pared to the standard factorization. Since the singularity associated with r 2 = is not 
even suppressed by the pion wave function, Sudakov effects are more crucial in Compton 
scattering [|T7| than in the case of form factors. The numerical outcomes of eq. ( P^ ) have 
been obtained and will be shown in section 5. 

Following the similar procedures, we derive the first invariant amplitude Hi. The 
extraction of Hi can be performed by contracting the two photon vertices with e^e' 1 ' 1 '. 
The derivation of Hi is much simpler than that of H 2 , because is orthogonal to all of 
the momenta Pi and g,. The contribution to the hard scattering associated with Hi from 
each diagram in fig. 3 is also listed in table. 1. The modified perturbative expression for 
Hi is given by 



Wi 



max 



(33) 



i=i 



Hi(s,t) = ^2 I dxidx 2 <f)(xi)<f)(x 2 ) I bdbT^(xi,s,t,b) 



x exp[-S(xi, 6, Q)} 



(34) 
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where 



f(i) = 8nC F (el + el)at a (wi) 
Hl (l-x 1 )(l-x 2 ) 

from the classes of fig. 3a-3c, and 

i(i) 



u s 

— I h4-2x-2y 

s u 



r\\b 



H 2 



l6irC F e u e d a s (w2) 
1 



(H -r,)K n ( V '|/-,|/; ) + e(r 2 )^H^ 



x 



+ 



1 



t 



Xi(l — Xi) X2{1 — X2) ^2(1 — XijS Xi{l — X2)U 



(35) 



(36) 



from the classes of fig. 3d-3e. The expressions for the Sudakov exponents and for wi are 
the same as before. 



5 Numerical Results and the Cross Section of 77 — > 

Based on the sum rules and the modified perturbative expressions for Hi and H2 in the 
previous sections, we compute the magnitudes of the two helicities. Sum rule predictions 



are obtained from eq. (15) with the substitution of sq = 0.6 and M 2 = 4 GeV 2 . The 



modified PQCD formula for H\ in eq. (0) is evaluated numerically, and H 2 is derived 
by H 2 = H — Hi, where the values of H have been given in || p| . Results of Hi 
and H 2 at different photon scattering angles 9* are shown in fig. 4a and 4b respectively, 
in which denotes the magnitude of Hi. Note that sum rule predictions for Hi are 
negative, and those for H 2 are positive. It is observed that sum rule results decrease more 
rapidly with momentum transfer \t\ [[Tj|, and have weaker angular dependence compared 
to PQCD ones. The PQCD predictions are always larger than those from sum rules at 
6* = 50° (—t/s = 0.6), and are always smaller at 6* = 15° (—t/s = 0.2) in the range 
4 < \t\ < 16 GeV 2 . It implies that large-angle Compton scattering might be dominated 
by perturbative dynamics. These two approaches overlap at 9* = 40° (—t/s = 0.5) and at 
\t\ = 4 GeV 2 , showing the transition of pion Compton scattering to PQCD. The transition 
scale is higher at smaller angles. Basically, the behaviors of Hi and H 2 are similar to that 
of H as shown in 0. 

With the knowledge of Hi and H 2 , we compute the cross section of pion Compton 
scattering. The expression for the differential cross section of pion Compton scattering in 
the Breit frame is derived in Appendix C: 

do" _ \Hi\ 2 + \H 2 \ 2 f s-u \ 3 
dcos#* 256vrt V s ) ' 1 ' 
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which yields the results exhibited in fig. 5. In the angular range we are investigating, the 
PQCD and sum rule methods predict opposite dependence on 9*: PQCD results increase, 
while sum rule results decrease, with the photon scattering angle. The reason for the 
difference is that the increase of the amplitudes with 9* is not sufficient to overcome the 
increase of the incident flux (see Appendix C). At a fixed angle, the differential cross 
section, similar to Hi and H 2 , drops as \t\ grows. Again, the transition scale is around 4 
GeV 2 for 9* = 40°. 

The differential cross section of pion Compton scattering from a polarized photon has 
been analyzed based on the standard factorization formula in [f[9|j . It is worthwhile to 
compare their predictions with ours from the modified perturbative formalism. Note that 
in [|19| the coupling constant a s is regarded as a phenomenological parameter and set to 
0.3, while we consider the running of a s due to the inclusion of radiative corrections, and its 
cutoff is determined by Sudakov suppression. Furthermore, our perturbative calculation 
is self-consistent in the sense that short-distance (small-6) contributions dominate. 

We concentrate on the two specific processes: 7r7t — > 7r7t and 7lvt — > 7r7t, where 7r 
(7l) denotes a photon with right-handed (left-handed) polarization. In our approach the 
amplitude of 7r7t — > 7 R 7r is derived in Appendix C as 



and that of 7lvt — > 7r7t is 



. . H\ — Ho 

M RR = - , (3* 



M LR = -*L + * . (39) 



The analysis in |19j was performed in the center-of-mass frame, in which the differential 
cross section is written as 

dcosfc' Szrrs 

with 9 the center-of-mass scattering angle, cos 9 = (t — u)/s. Substituting eqs. fl38|) and 
(|39"D into (|40|), we derive the modified PQCD predictions, evaluated at \t\ = 4 GeV 2 , 
whose dependence on 9 is shown in fig. 6a, along with the corresponding results obtained 
in |19[ using the Chernyak and Zhitnitsky wave function. It is found that the behavior 
of the differential cross section with respect to \t\ and 9 is similar to that in the Breit 
frame. The phase angles associated with A^rr and A^lr are shown in fig. 6b. Good 
agreement between the two predictions is observed in both of the figures, which justifies 
our perturbative calculation. Note that A^lr is real in the standard PQCD approach 
[|i~9"l , and its phase angle vanishes. 

The expressions presented in this paper can be further applied to another process, two- 



pion photoproduction 77 — > n + n~ . By a simple interchange of s and t in eqs. (|l^) and 
(|34T), we write down the sum rules and the modified perturbative expressions describing 
the process at intermediate Mandelstam variables. The motivation to study this process 



is that experimental data [1141 are available for the total invariant mass of the two pions, 
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M(7r + 7r~) = y/s, from 0.3 upto 1.5 GeV and for | cos^l < 0.6, 9 being the center-of-mass 
scattering angle. It is easy to check that the resulting sum rules for Hi involved in pion 
photoproduction are symmetric in 9 and tt — 9, as expected. 

However, restriction exists in the application due to the strong resonant contribution of 
/2(1270) appearing above s = 1 GeV 2 , which is not considered in our analysis. The region, 
in which it is possible to compare our predictions with data, is then limited to below the 
/2(1270) resonance. At such an energy scale, the PQCD formalism, even the modified one, 
can barely be reliable. Therefore, we will concentrate on the sum rule method. Another 
reason for favoring the sum rule formalism is due to the angular range of the data. It 
suffices to study the behavior of if; for < cos^ < 0.6, or, 0.2 < — t/s < 0.5, because 
of the angular symmetry mentioned above. In this range sum rule predictions are found 
to be dominant from the study of pion Compton scattering. On the other hand, the 
sum rule method does not work at very low energy, because the OPE of the correlator is 
not applicable when relevant scales go below 1 GeV 2 . Due to these constraints, we will 
investigate only the comparision of sum rule results with the experimental data in the 
region of M(tt + tt~) around 1 GeV and of | cos#| < 0.6. An agreement of our predictions 
with the data will justify the sum rule formalism given above. 

At s = 1 GeV 2 the asymptotic expressions for the perturbative spectral densities and 
the gluonic power corrections shown in section 3 are not appropriate, since s, t and the 
virtualities Si are of the same order of magnitude. It is easy to find that every term in 
the full series of pf ert , shown in Appendix B, are equally important at this scale. This 
fact requires us to employ the full expressions for pf ert and pf luon in the sum rule analysis 
of the process. Note that the quark power corrections (7? uark are exact. However, pf luon 
is too complicated to obtain its complete formula. The difficulty can be overcomed, if 
there indeed exists a largest stable region for Hi at some best so- Since the gluonic power 
correction should be about of the same order as the quark one, the best sq, obtained from 
the stability analysis based only on pf ert and Cj quark , is close to the exact value. As an 
approximation, we perform the stability analysis, say, of H 2 , by considering P2 Srt {s «-> t) 
given in Appendix B and C*2 Uark (s <-> t) at M{ti + ix~) = 1 GeV. The best Sq is found to 
be 0.3 GeV 2 , for which H2 does not vary much as M 2 > 4 GeV 2 for different 9. Once 
the approximate best so is determined, we compute Hi and if 2 in the large M 2 limit, 
the region where power corrections are negligible. By this means, the difficulty from 
the gluonic power corrections is avoided, and we need to evaluate only the perturbative 
spectral densities. Following the above procedures, we obtain Hi « —0.2 and H 2 ~ 0.1 
for different 9. We emphasize that at such a low energy scale power corrections must play 
an important role, and our results should be regarded as a rough estimation at most. 

Substituting the approximate values of Hi and H 2 into eq. fl40D , the total cross section 
(7(77 — > 7r + 7r~) at M(tt + it~) = 1 GeV for | cos#| < 0.6 is simply derived as 113 nb. The 
cross section at other values of M(tc + ti^) around 1 GeV can be computed in a similar way, 
and is found to be almost constant. Results are shown in fig. 7, along with part of the 
experimental data obtained by the Markll collobaration WM. It is obvious that our sum 



16 



rule estimation coincides with the data very well at M(tt + tt^) = 1 GeV. Above 1 GeV 
the data points rise rapidly due to the /2(1270) resonance as mentioned before. Below 1 
GeV the data show a slow falloff, and deviate away from our predictions, indicating that 
the sum rule method is not applicable at very low energy. Our formalism can be easily 
generalized to study another similar process 77 — > 7r°7r°, for which experimental data are 
also available 



6 Conclusions 

In this work we have extended the sum rule analysis to pion Compton scattering, compared 
their predictions with those from the modified factorization theorems, and shown that 
there is a clear overlap between the two approaches. We have given two individual sum 
rules and the modified perturbative expressions for the invariant helicity amplitudes to 
lowest order in a s . A detailed numerical analysis shows that the transition to PQCD in 
pion Compton scattering appears at \t\ = 4 GeV 2 and at 6* = 40° (in the Breit frame). 
This suggests that the sum rule and PQCD methods are complementary tools in the 
description of exclusive reactions, and can help locate their transition region by studying 
the power-law falloff of the corresponding amplitudes. Note that the stability of the sum 
rules is observed only when the modified phenomenological model is adopted. The values 
of the duality interval and the Borel mass are fixed at 0.6 and 4 GeV 2 respectively, values 
which are comparable to those commonly used in the sum rule analysis of pion form 
factor. 

We have also compared the modified PQCD predictions with those obtained by the 
standard perturbative calculation, in the center-of-mass frame ||19|| , and good agreement 



is observed. Note that our modified formalism does not involve a free parameter, and is 
relatively self-consistent. We have also shown, along the way, how to relate the helicities 
Hi to the description of pion Compton scattering from a polarized photon. As a last step, 
we have extended the sum rule formalism to the crossed process 77 — > 7r + 7r~ by simple 
crossing of s and t, and our predictions match the data of pion photoproduction. A more 
convincing justification of the formalism presented in this paper can be found from the 
study of other similar processes, such as proton Compton scattering [^T|. In our analysis, 
which is restricted to lowest order in a s , the sum rules for the helicities are real, and 
therefore, issues related to the perturbative and non-perturbative nature of the phases 
of Compton scattering cannot be addressed. We leave the discussion of these issues to 
future work. 
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Appendix A. Kinematics. 

This appendix is intended to provide a more detailed information on the methods 
employed in the calculation of the spectral densities of Compton scattering. Due to the 
great complexity involved in the analysis, we do not show every step of the computation, 
but illustrate the evaluation of some important integrals, especially those related to the 
new projection which lead to the sum rule for the invariant amplitude H\. Although the 
projector is a complex 4-vector, it is possible to prove that all the complex contributions 
cancel. The methods given in this appendix can be applied to a large class of Compton- 
type processes. The notations here are the same as in M, [TTJ. 

In the Breit frame of the incoming pion we have the Mandelstam invariants 

t= (P2-P1) 2 = s 1 + s 2 -2Q 2 Si> ' 2 



2Q 2 

glf2 

2Q 2 



u=( P 2-qi) 2 = 2Q 2 -s + ^ , (A.l) 



where Q 2 acts as a large parameter in the scattering process 

1 , 1 



Q 2 = ^(si + s 2 -t + 8) = ^(s + u + 6) , (A.2) 



with 

AQ A _^_ Sl s 2 
2Q 

Note that t = -2Q 2 at Sl = s 2 = 0. 



s 1 + s 2 -t) 2 -As lS2 = ^ o ^ 2 . (A.3) 



It is convenient to introduce light-cone variables for the momenta. We define 

qi = q+v + q^v + qi±, 

Pl = Qv + ^Lv, p 2 = + Qv, (A.4) 



where 

v = ^=(l,l,0 ± ), v= ^=(1,-1, Ojl) (A.5) 

satisfy the relations v 2 = v 2 = and v ■ v — 1. In our notation 

p± = i=(p ±p 3 ). (A.6) 
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The covariant expressions for qf can be easily obtained as 

{s-2Q 2 )(2Q 2 -s 2 ) 



Qi = 



2Q5 



_ (2Q 2 - Sl )(2Q 2 s- Sl s 2 ) 

The expressions for q^ are similar. 

In the frame specified above, all the transverse momenta are carried by the two pho- 
tons, which are on shell (qf, q\ = 0), and carry physical polarizations. The polarization 
vectors are defined by 

e (l)A = iV e (2)A r 



a „ A 



N x = e^Pp^R^ P a = Vip x + V2p x + _ R 

Vl=Pl-P2~S 2 , V2=Pl-P2-Su 

R = qi + q 2 , r = q 2 -q\, (A.8) 

with 

a = {2 /t) {y lPl ■ R + v 2 p 2 ■ R) . ( A.9) 

They satisfty the normalization conditions 

e W. gi = e w -g 2 = 0, 
e (0 . e 0) = . (A.10) 

These relations hold for all positive si and s 2 , whether or not they are equal, and for v 1 
and v 2 chosen as above. 

In this frame the projector n M = (n + ,n~ ,n±) is given by 



+ + ^fai* + g 2 ± )) /v 7 ^, 



n± = ^/=pf(9i + 92)± + ie (1) , (A.ll) 

which is a complex vector. 

All the integrals involved in the computation of the leading spectral functions corre- 
sponding to fig. la appear in the forms 

me,^,.)\ = /w,*-i,...o w ^-* ) ' ) ^- fc) ' ) . 

J {pi-k + qi) 1 

I'[f{k\ k ■ pi, ...)] = | rf 4 A;/(A; 2 , k ■ Pl , ...)5 + (k 2 )5 + ((p 1 - k) 2 )5 + ((p 2 - k) 2 ) , 

(A.12) 
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where / is a functin of arbitrary products containing the internal momentum k. The 
(^-functions come from the cutting rules applied to the quark lines except the top one. 
The components of k are then fixed at the values 

1+ = Q (2Q 2 s 2 - gigg) 
4Q 4 - Sl s 2 ' 

and k\ = 2k + k~. Hence, the only nontrivial integral needed to perform is the one 
associated with the polar angle 9 of k ± , which can be cast into the general form 

f 2n ,.sin ni 9 cos" 2 9 /A 

Tm,n 2 = d9— — . A.14 

v ' ; Jo A + Bcos9 K J 

The expressions for A and B depend on the diagrams we are considering. 

As an example we consider the following integral appearing in the evaluation of A\ ert : 



J[( e M • kfn ■ k] = I[(e ( l ] -k ± ) 2 (n+k- +n~k + )} 

a 



+/[(e? ) • k ± ) 2 (-^=q ± ■ k ± + te^ ■ k ± )} 



= \{n + k- + n-k + )klT(2, 0) - |gjf^T(l, 1) j 

xJ(Q 2 ,s u s 2 ) , (A.15) 
where J(Q 2 , si, s 2 ) is a jacobean given by 

As we have already discussed before, it is possible to show that the most significant 
contribution to the spectral functions comes from the region in which s\ and s 2 are close 
(si « s 2 ). 

Another useful integral for gluonic power correction corresponding to the diagram with 
a gluon from the vacuum attached to each of the lateral quark lines is 

i[f, mi , m2] . /^/^Mfa-^-^Mfe-*)'-^) (A17) 

J {Pi-k + qiY 

The relevant quantities are then obtained through the derivatives of this integral with 
respect to the two mass parameters mi and m 2 evaluated at mi = m 2 = 0. Using a 
coincise notation we define 

d d 

Id[S] = ^ a 7 [/> m l' m 2]Ui,m 2 =0 • (A.18) 

om 2 ami 
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Eq. ( |A. 17| ) can be computed in a similar way with the components of k fixed at 

~, + _ Q {2niiQ 2 - 2m 2 Q 2 - niis 2 + 2Q 2 s 2 - sis 2 ) 

" 4Q 4 - Sl s 2 ' 

q (2miQ 2 - mis x + m 2 s x + 2Q 2 s 1 - s 1 s 2 ) 1f1 v 

fc = 77^1 • (A. 19) 

4Q 4 - sis 2 

The angular integrals involved in this integral are more complicated, because A and B 
are given by irrational functions of Q 2 , s\ and s 2 , which, however, become rational when 
the masses are set to zero. 

Those integrals containing the factor k ■ qi in the numerator, either the type of 
eq. ([A. 12 ) or (|A.17 ), can be simplified by the following replacement: 



I[f, k- gi ]= Pl - gi I[f] + \ (/[/, (ft - k) 2 } - I'[f\) . (A.20) 

The expression for /'[/] has a strong similarity to that appearing in the calculation of 
form factors. For example, is the basic integral associated with the pion form factor: 

= J d*k5 + (k 2 )S + ((Pi - kf)5 + (( P2 - k) 2 ) = ^ . (A.21) 

In the /'-type of integrals the transverse momentum q±_ of the photon does not play any 
role as in the form factor case in the Breit frame. 

At last, for illustrative purposes we compute a simple example of eq. ( |A.17| ). Consider 



J[l, mi, 777,2] = J d 4 k 



S + (k 2 )S + ((p 2 - k) 2 - m 2 )6 + (( Pl - kf - m 2 ) 
(pi - k + <?i) 2 



_ 2 «) 3 1t(0,0) (A.22) 



with 



T(0,0) 



^ d6 



A + B cos 9 
2tt 



(A 2 - B 2 ) 1/2 ' 
A = s + m 1 -2k' + {p l +qi)~ -2k'~(p 1 + q x )' 



(A.23) 



B = yJSqt qi(2k'+k'- -mi) , 
A 2 - B 2 = (-4miQ 4 + 2m 2 Q 2 s + 4Q 4 s + 2m 1 Q 2 s 1 

—2m 2 Q 2 s\ — 2Q 2 ss x — 2m\Q 2 s 2 — 2Q 2 ss 2 

-misis 2 + ss!S 2 ) 2 / (4Q 4 - sis 2 ) 2 . (A.24) 
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The explicit expressions of most of the integrals are rather cumbersome and can not 
be listed here. The main feature of the calculation is the simplification provided by the 
use of the Breit frame, in which all the contributions are polynomials of Q 2 ,Si and s 2 . 
The spectral functions are symmetric in s\ and s 2 , at fixet s and t, as expected from 
time reversal invariance. For the sum rule investigation of Hi(s,t) at moderate energy 
scales, an expansion of the numerator of each term in the OPE up to leading virtualities 
of si and S2 is sufficient. We have examined this approximation numerically for different 
s and t, and found that the inclusion of terms with higher power of s\ and S2 gives only 
small modification. Based on the key integrals in eqs. (|A.12 ) and ( |A.17 ), it is possible 
to work out the expressions for the perturbative spectral functions, and the gluonic and 
quark power corrections given in sections 2 and 3. We refer in particular to || where the 
method of calculation of the power corrections has been developed. 



Appendix B. Perturbative Spectral Densities 

In this appendix we show the full expressions of the perturbative spectral densities 
for the invariant helicity amplitudes involved in pion Compton scattering. The derivation 
follows the points outlined in Appendix A, which is simplified drastically in the Breit 
frame of the pion. We summarize the results for p\ ert as 

pert _ _5 Rl(Q 2 , S, Si, Sg) , R , \ 

Pl ~ 24tt 2 Q 4 (2Q 2 - S ) 2 (4Q 4 - sis 2 f{2Q 2 s - Sl s 2 ) 2 ' 1 ' ' 

where 

15 

Ri{Q 2 , s, si, s 2 ) = an{s, s h s 2 )Q 2n . (B.2) 

n=0 

with 

9 9 

a = si s 2 

CLi = -2si 8 s 2 8 (4s + S! + s 2 ) 
a 2 = Assi 7 s 2 7 (7s + 3si + 3s 2 ) 

a 3 = -4si 6 s 2 6 (l2s 3 + 7s 2 si + 7s 2 s 2 - 8ssiS 2 - 4si 2 s 2 - 4siS 2 2 ) 
a 4 = 16si 5 s 2 5 (3s 4 + 3s 3 si + 3s 3 s 2 - 16s 2 sis 2 

— 8ssi 2 s 2 — 8ssiS 2 2 — 6si 2 s 2 2 ) 
a 5 = — 16si 4 s 2 4 (l3s 4 Si + 13s 4 s 2 — 20s 3 SiS 2 — 16s 2 Si 2 s 2 — 16s 2 SiS 2 2 

-40ssi 2 s 2 2 - 3si 3 s 2 2 - 3si 2 s 2 3 ) 
a 6 = 64si 3 s 2 3 (6s 5 si + 6s 5 s 2 + 22s 4 SiS 2 + 19s 3 Si 2 s 2 + 19s 3 s ; s 2 2 



22 



+ 19s 2 si 2 s 2 2 + 5ssi 3 s 2 2 + 5ss 1 2 s 2 3 + 2si 3 s 2 3 ) 
a 7 = -64 Sl 2 s 2 2 (4s 6 Sl + 4s 6 s 2 + 72s 5 sis 2 + 76s 4 s 1 2 s 2 + 76s 4 s 1 s 2 2 + 264s 3 siV 



+ 103s 2 si 3 s 2 2 + 103s 2 si 2 s 2 3 + 120ssi 3 s 2 3 + 12 Sl 4 s 2 3 + 12si 3 s 2 4 



a 8 = 



a 9 



256si 2 s 2 2 (l2s 6 + 30s 5 si + 30s 5 s 2 + 164s 4 SiS 2 + 98s 3 si 2 s 2 + 98s 3 SjS 2 2 
+ 19Qs 2 s 1 2 s 2 2 + A2ss 1 3 s 2 2 + A2ss 1 2 s 2 3 + 23si 3 s 2 3 ) 
-512sis 2 (8s 6 S! + 8s 6 s 2 + 72s 5 s 1 s 2 + 91s A s 1 2 s 2 + 91s A s 1 s 2 2 
+236s 3 si 2 s 2 2 + 92s 2 si 3 s 2 2 + 92s 2 si 2 s 2 3 + 84ssi 3 s 2 3 + 8si 4 s 2 3 + 8si 3 s 2 4 ) 
aio = 1024sis 2 (l2s 6 + 30s 5 si + 30s 5 s 2 + 134s 4 sis 2 

+98s 3 si 2 s 2 + 98s 3 sis 2 2 + 121s 2 si 2 s 2 2 + 27ssi 3 s 2 2 + 27ssi 2 s 2 3 + 2si 3 s 2 3 ) 
-1024s (4s 5 S! + 4s 5 s 2 + 72s 4 SiS 2 + 76s 3 Si 2 s 2 + 76s 3 SiS 2 2 
+204s 2 si 2 s 2 2 + 73ssi 3 s 2 2 + 73ssi 2 s 2 3 + 32si 3 s 2 ; 



a n = 



a 12 = 4096 (6s 5 si + 6s 5 s 2 + 37s 4 sis 2 + 19s 3 si 2 s 2 + 19s 3 sis 2 2 
+32s 2 Si 2 s 2 2 + 2ssi 3 s 2 2 + 2ssi 2 s 2 3 + 2sx 3 s 2 3 ) 

a 13 = -4096 (l3s 4 si + 13s 4 s 2 + 28s 3 sis 2 + 8s 2 si 2 s 2 + 8s 2 sis 2 2 
+8s Sl 2 s 2 2 + si 3 s 2 2 + si 2 s 2 3 ) 

ai4 = 16384s (3s 2 Si + 3s 2 s 2 + ssis 2 + Si 2 s 2 + Sis 2 2 ) 

ai 5 = -16384s 2 (si + s 2 ) . 

Obviously, the numerator Ri is a polynomial of the variables Q 2 , s and Sj. 
Similarly, the perturbative spectral density p 2 ert is given by 



pert _ _Jj i? 2 (Q 2 , S, Si, S 2 ) 

' - " 24tt 2 Q 4 (2Q 2 - s) 2 (4Q 4 - s lS2 f{2QH - Sl s 2 ) 2 ' 



where 



15 

R2{Q 2 , s, Sl , s 2 ) = K(s, si, s 2 )Q 2n (B.4) 

n=0 



with 



b = -si 9 s 2 9 
6i = 2si 8 s 2 8 (5s + si + s 2 ) 
h = -4si 7 s 2 7 (lOs 2 + 4ssi + 4ss 2 + sis 2 ) 
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6 3 = 4si 6 s 2 6 (20s 3 + lls 2 si + lls 2 s 2 - 6ssis 2 - 2si 2 s 2 - 2sis 2 2 ) 

6 4 = -16si 5 s 2 5 (5s 4 + 9s 3 si + 9s 3 s 2 - 20s 2 sis 2 - 10ssi 2 s 2 - 10ssis 2 2 - 8si 2 s 2 2 ) 
b 5 = 16si 4 s 2 4 (39s 4 si + 39s 4 s 2 + 4s 3 sis 2 -2s 2 si 2 s 2 -2s 2 sis 2 2 

-62ss! 2 s 2 2 - lls! 3 s 2 2 - lls! 2 s 2 3 ) 
b 6 = -64si 3 s 2 3 (l8s 5 si + 18s 5 s 2 + 86s 4 siS 2 + 65s 3 Si 2 s 2 + 65s 3 S!S 2 2 

+78s 2 si 2 s 2 2 + 16ss! 3 s 2 2 + 16ss! 2 s 2 3 - S! 3 S 2 3 ) 
b 7 = 64si 2 s 2 2 (l2s 6 si + 12s 6 s 2 + 216s 5 s 1 s 2 + 236s 4 si 2 s 2 + 236s 4 sis 2 2 

+760s 3 si 2 s 2 2 + 273s 2 si 3 s 2 2 + 273s 2 si 2 s 2 3 + 290ssi 3 s 2 3 + 32s! 4 s 2 3 + 32si 3 s 2 4 ) 
b 8 = -256si 2 s 2 2 (36s 6 + 90s 5 Si + 90s 5 s 2 + 464s 4 S!S 2 + 286s 3 si 2 s 2 + 286s 3 Sis 2 2 

+516s 2 S! 2 s 2 2 + 106ssi 3 s 2 2 + 106ssx 2 s 2 3 + 59s! 3 s 2 3 ) 
69 = 512sis 2 (24s 6 si + 24s 6 s 2 + 216s 5 sis 2 + 265s 4 si 2 s 2 + 265s 4 sis 2 2 

+676s 3 si 2 s 2 2 + 266s 2 si 3 s 2 2 + 266s 2 si 2 s 2 3 + 231ssi 3 s 2 3 + 20si 4 s 2 3 + 20si 3 s 2 4 ) 

610 = -1024sis 2 (36s 6 + 90s 5 S! + 90s 5 s 2 + 414s 4 SiS 2 + 286s 3 Si 2 s 2 + 286s 3 SiS 2 2 

+386s 2 S! 2 s 2 2 + 86ssi 3 s 2 2 + 86ssi 2 s 2 3 + 13si 3 s 2 3 ) 

611 = 1024 (l2s 6 si + 12s 6 s 2 + 216s 5 sis 2 + 236s 4 si 2 s 2 + 236s 4 sis 2 2 

+660s 3 S! 2 s 2 2 + 233s 2 S! 3 s 2 2 + 233s 2 Si 2 s 2 3 + 150ss! 3 s 2 3 + 10si 4 s 2 3 + 10si 3 s 2 4 ) 
b 12 = -4096 (l8s 5 si + 18s 5 s 2 + llls 4 sis 2 + 65s 3 si 2 s 2 + 65s 3 s 1 s 2 2 

+ 108s 2 si 2 s 2 2 + 12ssi 3 s 2 2 + 12ssi 2 s 2 3 + 8si 3 s 2 3 ) 
613 = 4096 (39s 4 si + 39s 4 s 2 + 84s 3 SiS 2 + 30s 2 S! 2 s 2 + 30s 2 Sis 2 2 

+22ssi 2 s 2 2 + si 3 s 2 2 + si 2 s 2 3 ) 

b 14 = -16384 (9s 3 si + 9s 3 s 2 + 2s 2 sis 2 + 2ssi 2 s 2 + 2ssis 2 2 - si 2 s 2 2 ) 
615 = 16384s (3ssi + 3ss 2 — 2sis 2 ) . 

Based on the above formulas, it is easy to confirm that those terms with higher power of 
Si and s 2 are negligible, and the asympotic expressions in eq. (|T5|) are justified. The full 
expressions for pf CTt listed here are useful in the analysis of pion photoproduction. 

Appendix C. The Scattering Amplitude and Cross Section 

In this appendix we derive the expressions for the amplitudes of pion Compton scatter- 
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ing from a polarized photon and the cross section appearing in section 5. The differential 
cross section of Compton scattering is defined by 



d(7 = \^l d Q , (C.l) 

where M. is the scattering amplitude, F is the incident flux, and dQ is the phase space 
of the final states, 

d Q = (2.)^>( P1 + * - „ - *) ( J;1 Ek(2 ^ • (C2) 

F is defined by 

F = \v Pl ~ v qi \2E Pl 2E qi =4p.q = 2s, (C.3) 

with v p = p/Ep the velocity of the incoming particle. It is easy to observe that F increases 
with the photon scattering angle. Combining eqs. (|C2 ) and ( |U.3| ) we obtain the general 
expression 

da= \^lL dcos (C.4) 



with 9 the center-of-mass scattering angle. The above formula is exactly eq. (f40|) in section 
5, which is employed in the comparision of our predictions with those from ref. [19] and 
with experimental data of 77 — > 7r + 7r~ [[14]]. 

Substituting cos 9 = (t — u)/s into eq. (|U.4|) , we have the expression which is invariant 
in both of the center-of-mass and Breit frames. Then eq. (|C.4|) can be easily converted 
into the one in the Breit frame using the relation sin 9*/2 = —t/(s — u) with 9* as defined 
before. We have 

dcr \M\ 2 f s — u\ 3 
dcos#* ~ 1287rt 
where the scattering amplitude M. is given by 



(C.5) 



M = M^ 1T e*» T (C.6) 

with defined by eq. (^), and €t the polarization vector of the photon in the state 
T. Inserting \A4\ 2 = {\Hi\ 2 + \H2\ 2 )/2 for pion Compton scattering from an unpolarized 
photon into the above formula, we obtain eq. (|37|) . 

In the case involving a polarized photon, the scattering amplitude is given by 

-Mrr = (C.7) 

for the process 7r7t — > 7r7t, and 

^lr = M^ L e- (C.8) 



25 



for 7l7t — > 7r7t. To compute AIrr and AIlr explicitly, we assign the following momenta 
in the center-of-mass frame ||19|| : 

Pl = P(l, 0,0,-1), qi = P(l, 0,0,1), 

p 2 = P (1, — sin 0, — cos^) , g 2 = -P(l, sin#, 0, cos#) , 

ei^r = -^(0,-l,±i,0), e 2Rii = -^=(O,-cos0,±i,sin0) . (C.9) 

A straightforward calculation leads to the relations 

e {1) • e 1R = -e« • e 1L = -e« • e\ R = ± , 

e (2) • e lR = e< 2 ) ■ e 1L = -e® ■ e* 2R = --L , (CIO) 
based on which it is then a simple matter to work out eqs. fl38|) and (j39|). 
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Table 1. The expressions of the hard scatterings T H and corresponding to the 
diagrams in fig. 3. Here we define 
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Figure Captions 

Fig. 1 Lowest-order diagrams for pion Compton scattering in the sum rule approach. 

Fig. 2 Variation of (a) Hi and (b) H 2 with M 2 for (1) s = 0.7, (2) s = 0.6, and (3) 
s = 0.5 GeV 2 . 

Fig. 3 Lowest-order diagrams for pion Compton scattering in the perturbative QCD 
approach. 

Fig. 4 Dependence of (a) and (b) |£||i?2| on \t\ derived from the modified PQCD 

(real lines) and from QCD sum rules (dashed lines) for (1) —t/s = 0.6 (8* = 50°) (2) 
-t/s = 0.5 (0* = 40°), and (3) -t/s = 0.2 (9* = 15°). 

Fig. 5 Dependence of da/dcos6* on \t\ derived from the modified PQCD (real lines) and 
from QCD sum rules (dashed lines) for (1) —t/s = 0.6 (2) —t/s = 0.5, and (3) —t/s = 0.2. 
Note that the curve (1) from sum rules is shown by a long dashed line. 

Fig. 6 Dependence of (a) S 3 da/dcos6, and of (b) the corresponding phase in degrees 
on cos8 from the modified PQCD (real lines) and from ref. |19| (dashed lines) for (1) 
7r7t — ► 7r7t and (2) 7l7t — ► 7r7t. Note that S = sin(6 l /2) here. 

Fig. 7 Dependence of 17(77 ~~ > 7r + 7r~) on M(tt + 7t~) derived from QCD sum rules. Part 
of experimental data adopted from |T^] are also shown. 
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